Effect of conformations on charge transport in a thin elastic tube 
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We study the effect of conformations on charge transport in a thin elastic tube. Using the Kirchhoff 
model for a tube with any given Poisson ratio, cross-sectional shape and intrinsic twist, we obtain a 
class of exact solutions for its conformation. The tube's torsion is found in terms of its intrinsic twist 
and its Poisson ratio, while its curvature satisfies a nonlinear differential equation which supports 
exact periodic solutions in the form of Jacobi elliptic functions, which we call conformon lattice 
solutions. These solutions typically describe conformations with loops. Each solution induces a 
corresponding quantum effective periodic potential in the Schrodinger equation for an electron in 
the tube. The wave function describes the delocalization of the electron along the central axis of 
the tube. We discuss some possible applications of this novel mechanism of charge transport. 
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I. INTRODUCTION 

Nanotubes and nanowires have attracted considerable 
attention due to their their potential technological appli- 
cations [1]. Carbon nanotubes and DNA are well known 
examples. One important problem is to understand their 
mechanical properties. Another is to gain some insight 
into the effect of quantum confinement on electronic 
transport in such nanostructures. 

One striking observation concerning the DNA molecule 
is that its central axis can take on various curved con- 
formations. Similarly, nanotubes need not be always 
straight. Helix-shaped nanotubes have been observed by 
Volodin et al^. One way to theoretically model these 
twisted thin tubes (wires) would be to assume that they 
are elastic filaments described within the framework of 
the Kirchhoff rod model [3], derive the various possible 
curved conformations that can arise from it, and see 
if they are in accordance with the shapes actually 
observed. This would also help in studying the elastic 
properties of these structures within this model. Q The 
natural question that arises is how the confinement of a 
particle to a curved path would affect its transport. 

Now, the problem of quantum transport of a free 
particle in a thin curved tube has been studied by several 
authors [1, d, 0|- They have shown that the curved 
geometry of the tube essentially induces a potential in 
the path of the particle, thereby affecting its transport 
properties. 

In a recent paper Q, we had analyzed the case when 
the thin tube is made of elastic material, and applied it 
to electron transport in a biopolymer. We analyzed the 
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statics and dynamics of the thin elastic tube using the 
Kirchhoff model f^. Under certain conditions, the curva- 
ture function of the tube supports a spatially localized 
traveling wave solution called a Kovalevskaya solitary 
wave d, • We showed that this solution corresponds 
to a conformation which is a localized twisted loop 
traveling along the tube Q. The localized bend induces 
a quantum potential well in the Schrodinger equation 
for an electron in the tube [1, @, 0] j which traps it in the 
twisted loop , resulting in its efficient motion, without 
change of form, along the polymer. Our result formalizes 
the concept of a conformon that has been hypothesized 
in biology. 



The motivation for the present paper is as follows. 
Firstly, the analysis of the Kirchhoff model given in [^IToj 
had assumed that the elastic tube had no intrinsic twist. 
It was therefore not strictly applicable to either to DNA 
(which has an intrinsic twist of about lOdeg/A, even 
in the straight relaxed state) or to the class of intrinsi- 
cally twisted nanotubes. In what follows, we present a 
more general analysis of the Kirchhoff model, by incor- 
porating an intrinsic twist in the tube, and point out the 
nontrivial role that it plays in determining the geometri- 
cal shape of the axis of the tube. Further, our analysis 
is valid for any general value of the Poisson ratio. Sec- 
ondly, we show that the nonlinear differential equation 
for the curvature function obtained from the static Kirch- 
hoff equations can, in addition to the localized solution 
discussed in also support a class of spatially periodic 
solutions in the form of Jacobi elliptic functions. We call 
them conformon lattice solutions. Each such solution in- 
duces a quantum periodic potential. Finally, we show ex- 
plicitly that the Schrodinger equation supports an exact 
solution, which corresponds to the electron getting delo- 
calized along the axis of the quantum wire, even in the 
static case. Thus this mechanism for electron transport 
in a quantum nanotube is distinct from the conformon 
mechanism discussed in Ref. B. 



2 



II. THE KIRCHHOFF MODEL 

We consider a thin elastic rod [3] whose central axis is 
described by a space curve R(s,i). Here, s denotes the 
arc length of the curve and t the time. Let the rod remain 
inextensible with time evolution, with the unit tangent 
to its axial curve given by t = R^. (The subscript s 
denotes the derivative with respect to s). In the plane 
perpendicular to t, we define as usual , two orthogonal 
unit vectors n and b, where the principal normal n is the 
unit vector along t^, and the binormal b = (t x n). The 
rotation of the Frenet frame (t, n, b) as it moves along the 
curve, is given by the well known Frenet-Serret equations, 
ts = kn, Us = —kt + rb and hg = — Tn, where the 
curvature k = \ts\ and the torsion t — t ■ {tg x tgs) / k'^ . 

For elastic tubes, instead of n and b, it is more natural 
to use 'material' orthogonal unit vectors di and d2 lying 
on the cross-section of the tube, perpendicular to its axial 
curve. For e.g., they could lie along the principal axes 
of inertia of the cross-section. Let di make an angle (f> 
with n. Hence we get di = ncoa cf) + bsin0 and d2 = 
—nsincj) + bcos0. Denoting the tangent t by da, and 
using the Frenet-Serret equations, it is easy to show that 
the rotation of the 'material' frame (di, d2, d^) along the 
curve is given by 

di s = k X d,, (1) 
where i = 1, 2, 3, and the vector k is given by 

k = fcidi fczdz -f fcgdg. (2) 

Here 



In the equations above, the parameter a is a measure of 
the bending rigidity of the cross section of the rod. It 
is the ratio /1//2 of the two moments of inertia of the 
cross-sections of the nanotube in the directions di and 
d2 . These are conventionally oriented such that Ii < l2- 
Hence < a < 1. (For a circular cross section, a — 1.) b 
is a measure of the twisting rigidity of the rod. It is given 
in terms of a and the Poisson ratio a, which is a measure 
of the change in volume of the rod as it is stretched: 

6 = 2a/[(l + a)(l + a)]. (8) 

In Refs. [1, [l^l, the term involving fcg in Eq. ([7]) is 
absent. Hence the analysis given there is valid for tubes 
without an intrinsic twist. In the next section, we analyze 
the Kirchhoff' equations (g]) and © with fc!^ ^ 0. Our 
analysis will be valid for general values of a and a. 

III. ROLE OF THE INTRINSIC TWIST 

We first analyze the static conformations. Substituting 
Eqs. dH) and ^ in Eqs. (g]) and (O, we get 



9i,s + fe.g3 - ^352 = 0, (9) 

92,8 + ksgi - kigs = 0, (10) 

53, s + ^152 - ^251 = 0, (11) 

92 = ki^s + (6 - a)fc2fc3 - bk% , (12) 



(fci, fe, fcs) = (fcsin0, kcoscj), T + (f>s). (3) 

Let the internal elastic force (or tension) and the to- 
tal torque, that act on each cross section of the rod be 
given by g and m respectively. The Kirchhoff equations 
that result from the conservation of linear and angular 
momentum at every point s are given in dimensionless 
form by 

Ss = Rtt (4) 

and 

+ ds X g = adi x di^ft + d2 x d2,tt , (5) 

where the subcript t stands for the derivative with respect 
to time. The internal force is written in a general form 
as 

g = 5idi -I- .g2d2 + .g3d3 . (6) 

We take into account the intrinsic twist fcg of the elastic 
tube, by writing the constitutive relation for the internal 
torque m as 



= -ak2,s + {b- l)/ci/c3 - bk^ki , (13) 

and 

bk^.s + {a- I)fci/C2 = 0. (14) 

As is clear from Eq. fS]), the equations above represent a 
set of nonlinear coupled differential equations involving 
k, T and cf). 

To understand the role played by the intrinsic twist, 
we proceed as in Refs. 8, 10], and look for solutions that 
correspond to </) = ^nTr, n = 0, 1, . . . . Using this in Eq. 
(Ill), we see that Eq. yields 

T = To . (15) 

Thus the torsion of the tube, which is a measure of its 
nonplanarity, is an arbitrary constant, to be determined 
consistently. We focus on nontrivial conformations with 
a nonvanishing tq. Two cases arise in the analysis of 
Eqs. ©-(HSl): 

Case (i): <f) = jir, j = 0,1; ki — 0; k2 = {—yk. 
We find in this case 



m = fcidi -I- afc2d2 -I- ^(fcs - fc3)d3. 



(7) 



(6 - 2a) To^bk° 



(16) 
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Using Eq. ([8]) in this equation, we get 

To = -kl/[a + a{a + l)]. (17) 

The internal elastic force is found to be 

g=(-l)^a(rofcd2-A;,di) + (C-iafc2) , (18) 

where the constant C represents the tension in the rod 
when it is straight. 

Case (ii): = (j + ^)tt , j = 0, 1 ; fci = (-)Jfc ; fca = 0. 
We now obtain 

(6-2)ro = 6fc^ (19) 
As before, using Eq. ([H]) in this equation yields 

To = -afc5/[l + a(a + l)]. (20) 
The elastic force is determined as 



equation for the curvature k in Cases (i) and (ii). They 
are seen to have the same form 



g = (-l)-' [(rofcdi+fc, d2) + (C- 



iP) 



(21) 



We now show that the intrinsic twist indeed plays a 
nontrivial role in determining the torsion tq of the con- 
formation of the elastic tube. 

First, if /eg = 0, then Eqs. (HH) and ^ yield b ^ 2a 
and 6=2, respectively, because tq 7^ 0. Putting in 
these two values of b in Eq. ([8]) successively, we get 
a = — 1/(1 + a) and a = — a/(l -I- a), respectively. 
Since < a < 1, it follows that the Poisson ratio a 
has to be negative^ and in the range — 1 < cr < — i and 
— i < cr < 0, respectively, for the two cases. Although 
thermodynamic stability arguments [l^ merely restrict a 
to the range — 1 < cr < ^, and while a can indeed be 
negative for some biopolvmers [Toj . for most elastic me- 
dia one finds that < cr < i. 

Thus, for both cases (i) and (ii), setting — deter- 
mines the numerical value of cr, which turns out to be 
negative for any a. Further, tq can take on any arbitrary 
value in this instance. 

In contrast, if fcg 7^ 0, the corresponding torsion tq is 
not arbitrary, but is determined in terms of fcjj, a and 
cr (which are material properties), as we may expect on 
physical grounds. This dependence can be seen from Eqs. 
(fT7| and ([20| . These equations also show that positive as 
well as negative values of cr are allowed in this instance, 
which is a desirable feature. We may note that when 
CT > 0, the torsion tq and the intrinsic twist /cg have 
opposite signs. In addition, on imposing the condition 
< a < 1, both the equations (|17p and (PO)) lead to the 
same inequality, cr < — fcg/ro < 2cr + 1. 



IV. EXACT CONFORMATIONS 

In the last section, we found the solutions for the tor- 
sion Tq. Using Eq. ^ in Eq. dU, and Eq. ^ in Eq. 
(llOp , respectively, we can derive the nonlinear differential 



kss + \k^ = (C2-T2)fc, 



(22) 



where C2 = C /a for Case (i), while C2 = C for Case (ii). 

Note that Eq. (p2)) has the same form as that obtained 
earlier [1, [l^ , for the case fcg = as well, but with the 
important difference that tq is not arbitrary any more, 
but depends on k^ as shown in the last section. Eq. 
has a solution of the form 



k{s) = 2{C2 - To2)i/2 sech [(C2 - ro2) 



^l/2 



(23) 



for (C2 - r2) > 0. 

This represents a static conformon. Now, due to scale 
and galilean invariance Q of the dynamic Kirchhoff equa- 
tions (|4]) and ([5]) , traveling wave solution are possible for 
k, with s replaced by (s — vt) in Eq. psp . Here v rep- 
resents the velocity of the wave. This leads to a moving 
conformonQ, which is like a solitary wave. 

Now, in general, we find that the nonlinear differential 
equation (|^^ supports solutions of the form 



M^,'^) = 2,/§f^dn(y%^.,.), (24) 



2-k2 



where dn is the usual Jacobi elliptic function, with mod- 
ulus n in the range < k < 1 [16]. For k = 1, the 
solution ([M]) becomes ([25)1 . which was considered in Ref. 
8] . But in contrast to that solution, which was spatially 
localized, the solution Eq. ((24|) is a spatially periodic 
function, with a finite period 2K{k) for all k 7^ 1. Here 
K{k) is the complete Jacobi integral of the first kind, 
which tends to infinity for k — 1. Thus we call Eq. 
for all K 7^ 1, a conformon lattice solution. 

In Figs. 1 to 5, we have displayed the conformations of 
a thin elastic tube (rod) that correspond to the confor- 
mon lattice solutions for the curvature, given in Eq. (|24p . 
for various values of k, and constant torsion tq = 0.5. In 
all the plots, -y/ (C2 — to2) has been set equal to unity, 
since its inverse merely scales the length. 

For K = 1, the conformation is given in Fig. 1, and it 
corresponds to a single loop or conformon. Its shape is 
(as expected) different from that obtained in Ref. ( 8]) 
for To = 1. 

We find that even a slight decrease to k = 0.995 leads 
to a fairly large change in the conformation, with two 
loops, and a partially complete loop, as seen in Fig. 2. 
From K = 0.995 to k = 0.75, there are steady changes 
in the conformation. In contrast, below n — 0.75, the 
conformational changes are much slower as n varies, as 
seen by comparing Fig. 3 (k = 0.75) with Fig. 4 (k = 
0.25), which differ very slightly. 

For K = 0, the curvature becomes a constant indepen- 
dent of s. Hence this conformation represents a structure 
in which the axis of the tube is coiled into a helix, with a 
constant curvature and torsion. This is displayed in Fig. 



FIG. 1: Conformation corresponding to curvature k as given 
in Ea. (|24p for k = 1. Note the localized twisted loop formed 
by the axis of the tube. 




FIG. 2: Conformation corresponding to curvature k as given 
in Ea. (|24|) for k — 0.995. Note the appearance of two loops 
and a partially complete loop. 



FIG. 3: Conformation corresponding to curvature k as given 
in Ea. (l24p for k — 0.75. There are three twisted loops of 
unequal sizes, and a bend. 




FIG. 4: Conformation corresponding to curvature k as given 
in Ea. (|24|) for k — 0.25. There are three twisted loops of 
unequal sizes, and a bend. 



5. We parenthetically remark that while this helical tube 
was essentially the only conformation considered in ^4] in 
the context of the Kirchoff model, we see that a wide 
class of coiled conformations corresponding to nonzero k 
get supported by the model. 

As explained below the solution Eq. for the con- 
formon, it is possible to have traveling wave solutions 



(231) with s replaced by (s — vt) . 

We conclude this section with the remark that the solu- 
tions for the curvature (Eq. ((24)) ) will also be applicable 
for a closed tube of length L satisfying periodic boundary 

conditions. We get Ly 2-k° ~ 2m/f (k), where m is 
an integer. 
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FIG. 5: Conformation corresponding to curvature k as given 
in Eq. (|24|l . which becomes a constant for k = 0. There are 
four identical loops, representing four full windings of a helix. 



V. QUANTUM MECHANICAL CHARGE 
TRANSPORT IN A TUBE 

In this section, we consider the implications of the con- 
formon lattice solution P^ on electron transport along a 
nanotube. We recall some salient steps given in Q, where 
we discussed the conformon, for the sake of completeness. 

As has been shown @, 0|i a thin tube with curva- 
ture k(s) and torsion tq, induces a quantum potential 
Veff{s) = (?i^/2/x)[-fc2/4 + To2/2] for an electron mov- 
ing on it. Here fi is the mass of the electron. In the time- 
dependent Schrodinger equation with the above poten- 
tial, we can eliminate the constant torsion tq using a sim- 
ple gauge transformation ipi — -0(3, i) exp(— i?iTo2t/4/i). 
This leads to 



iTT'-Q-^i'{s,t) 



h2 / 82 

2^ vas2 



k2{s) 



i^{s,t) (25) 



Using the transformation (t, s) — > (^u, v^si), Eq. 
25l) becomes 



k2 



(26) 



where k = fc(si), and the subscripts si and u stand for 
the partial derivatives and 

Looking for stationary solutions of the Schrodinger 
equation (|26p in the form 



ip{si, u) ~ fc(si) exp{~iE u), 



(27) 



we get 



ksisi + 



fc3 



Comparing this with Eq. 
must have 



-E k. (28) 
for self-consistency, we 

(29) 



which is negative. We write Eq. ((28)l in the form of a 
time-independent Schrodinger equation with a potential 
term V 



V{s,K)k = Ek. 



(30) 



On using Eq. ([24ll . we see that the potential is given by 



F(s,k) = -2 



(g2-rg) 



dn2( 



C2 - ro2 
2-k2 



,k), (31) 



which is negative as well. Since <hi^{s,K) < 1 for all s 
and K, the minimum value of the potential in Eq. (|28p is 



V„. 



(C2-ro2)/[l 



(32) 



Thus we have the inequality > -B > Vmin- On using 
Eq. (IIH) and Eq. ((291) in Eq. (H?]) , the exact stationary 
solution of the time dependent Schrodinger equation for 
the electron (Eq. (|26p) is of the form 



V'(si, k,m) = 2\ 



'C2 



^dn(. 



IC2 



■f- si,k) expi (C2-Tq)u. 



(33) 



r 



Hence we see that the conformon lattice solution, a spa- 
tially periodic solution for k{si) given in Eq. (|24p . which 
determines the conformation of the elastic tube, is also 
the amplitude of the electron wave function in (|33p . 
Note that when n goes to 1, the amplitude of the above 



solution reduces to Eq. which corresponds to the 

static conformation in Fig. 5. This would in turn localize 
the electron at the center of the loop. 

However, as we explained below Eq. (|23p . a dynamical 
solution with s replaced by (s — vt) is possible for the 
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KirchhofF equations ^ and ([5]), which would describe a 
moving loop. Turning to electron transport on a tube, 
we showed in |8| that the quantum effective potential 
induced by the above (moving) curvature traps the elec- 
tron, which then travels with it along the tube. This 
is the conformon mechanism described in Ref. [8] and 
applied to a biopolymer.. 

Now, we see that the conformon lattice solution for the 
curvature (Eq. to a different type of charge trans- 

port mechanism. While for n = 1, it was the motion of 
the conformon that led to electron transport, we see that 
for K 7^ 1, even the static conformon lattice contributes 
to electron transport. This happens because the electron 
wave function Eq. corresponds to a spatially delocal- 
ized state, due to the finite periodicity of dn . Indeed, Eq. 
()33|) shows that since for any k, dn (s, k) does not vanish 
for any s, the probability of finding the electron on the 
polymer, dn^(s,K), is non- vanishing for all s. Note that 
the period of dn (s, k) increases as k increases, and tends 
to infinity as k tends to 1. Further, dn(s,K) itself goes 
to unity for all s as k tends to zero. The conformon lat- 
tice mechanism therefore contributes towards enhancing 
electron transport along the nanotube as k is decreased 
from 1 to 0. 



VI. DISCUSSION 

With regard to applications, let us first consider car- 
bon nanotubes. It has been shown using Green Function 
techniques in the context of an " arm chair" carbon nan- 
otube that bending of a nanotube increases its electrical 
resistance [TtI. This agrees with our finding that a single 
bend (as in Fig. 1) would cause localization of the elec- 
tron wave function. Thus although the phenomenological 
model we have presented is a highly idealized continuum 
model in the sense that it ignores details such as elec- 
tronic structure, it is able to capture the essence of the 
curved geometric effects quite succinctly. Its merits in- 
clude the use of a realistic elastic model that takes into 
account an intrinsic twist and a possible asymmetry in 
the bending rigidities. We have obtained a helical con- 
formation (Fig. 5) which has been observed experimen- 
tally. While the fact that a helical solution per se can 
be supported by the Kirchhoff model is well known p^ . 
here it emerges as a special case of an elliptic function 
solution. It would be interesting to experimentally mea- 
sure its electrical resistence as a function of the number 
of coils. We conjecture that the resistance of nanotubes 
where the coils appear very close to each other may be 
less than that in which the coils are far apart, since the 
delocalization of the electron will be very efficient in the 
former case. Hence design of suitable experiments to ver- 
ify the above would be of interest. 

More importantly, we have shown that in addition to 
the helical conformation, the KirchhofF model supports 
shapes such as those given in Figs. 1 to 4, which have 
constant torsion but varying curvature. While we have 



seen in [8[ that a static loop conformation (Fig. 1) lo- 
calizes the charge on the loop, the new multiloop solu- 
tions (Figs. 2 to 5) are shown to lead to delocalization 
of the wave function, which provides a novel mechanism 
for charge transport. 

Conformations corresponding to closed tubes are also 
possible. First, the simplest case of a circular ring (and 
hence planar) solution is easily obtained by setting the 
curvature /c to be a constant Kq, and the torsion tq = in 
the basic equations ([9]) to p4)) . We obtain = \J1Ci. 
This is the same as the solution of fc found from Eq. 
(1^ by inspection in this case. In addition, closed tube 
conformations for the general elliptic function solution 

([M]) exist when its length L is such that L\l is 
an integral multiple of its period , 2K{k), with periodic 
boundary conditions R(i) ~ R(0). 

In the context of DNA, Shi and Hearst pj| have anal- 
ysed the KirchhofF equations, by treating DNA as a thin 
elastic rod with a circular cross-section, to obtain exact 
solutions of curvature and torsion. They have plotted 
closed tube configurations of a supercoiled DNA in the 
form of a toroidal helix. Our analysis and plots are for 
open tubes with non-circular cross-section, the helix be- 
ing a special case. 

Over the years, physicists have used various types of 
experiments to measure DNA conductivity, but the re- 
sults contradict each other. 20] Bioscientists have em- 
ployed other types of methods[21] such as oxidative dam- 
age in DNA, to study charge transport. They have firmly 
established that there is definitely an efficient trans- 
port of electrons (or holes) through the DNA base pair 
stack, which proceeds over fairly long molecular distances 
(~ 200 A) under ordinary thermal conditions. Since the 
central axis of DNA can take on curved conformations, 
especially when it interacts with various proteins during 
replication and transcription, the quantum mechanical 
motion of the electron is not always along a straight line, 
but rather, on a curved path. This curvature perhaps 
plays a role in charge transport, since the sensitivity of 
transport to sequence-dependent conformations as well 
as local flexibilty and distortions of the soft DNA chain 
have been noted in experiments. (2ll. [2^ 

It is indeed not easy to measure coherent charge trans- 
port in DNA, because of effects due to temperature, dis- 
order, dissipation and interaction with molecular vibra- 
tions, etc |20|. In addition, issues connected with the 
source of the charge, whether it is transport between the 
leads or due to a donor-acceptor scenario, will need care- 
ful experimental scrutiny. But we believe that irrespec- 
tive of how the charge gets injected, given the same ambi- 
ent conditions like temperature, a coiled elastic nanotube 
induces a quantum potential in its path, while a straight 
tube of the same length does not, and this should lead to 
observable effects on charge transport. 

We note that the formation of DNA loops (such as 
in Figs. 1 to 5) has been frequently observed [23|, and 
is considered to be an important biological mechanism 
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for regulating gene expression. Thus the exact confor- 
mations that we have obtained using the clastic model 
indeed seem to have physical significance. In view of 
the above, systematic experiments should be designed to 
study the dependence of charge transport on DNA con- 
formations. 

We suggest that in such experiments, long molecules 
(with length greater than the persistence length of 50nm) 
which arc not stretched out or attached to a surface, but 
rather in their more natural, free curved conformation 
(with loops, if possible) should be used to investigate 
for what kinds of conformations, there would be either 
enhancement or depletion of charge migration. 

In summary, as a result of an intricate interplay be- 



tween classical elasticity of the nanotube and quantum 
motion of the electron in it, certain conformations with 
constant torsion arc shown to lead to a class of corre- 
sponding exact electronic wave functions. This suggests 
a charge transport mechanism that underscores the role 
played by the curved geometry of the axis of the tube, 
and in particular, the nonlinear differential equation sat- 
isfied by the curvature function. Hence curved geometry 
should be taken into account in order to obtain a full 
understanding of charge conduction in nanotubes. Our 
work should be regarded as a first step in this direction. 

RB thanks the Council of Scientific and Industrial Re- 
search, India, for financial support under the Emeritus 
Scientist Scheme. 
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